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ABSTRACT 


Remote sensing is an advanced technique that can be used to detect 
hydrocarbons on land, and its advantages include rapid and low-cost 
detection. The detection of hydrocarbons on the surface is essential for 
environmental monitoring and the purpose of exploration. For oil and gas 
companies, the ability to locate hydrocarbon spill sites is important 
information for the success of future exploration wells. This study aimed to 
determine the quantities of hydrocarbons mixed with the soil as an indicator 
of its accumulation in the subsurface. In this paper, temperature data 
obtained from the landsat-8 satellite were used to determine the extent of the 
oil spill in the area using spatial interpolation and gradient techniques. 
Kriging interpolation is based on statistical models (i.e. including 
auto-correlation) geostatistical techniques can produce a forecast surface and 
also provide accurate forecasts. The results showed that there are 60 sites of 
thermal anomalies and temperature values are lowest and highest 23.2 °C to 


91.11 °C, in the study area. Three different sites were chosen in terms of area 
and effect. They are located in the northern and southern Rumaila and 
Zubair fields. The location of the oil spills was determined with high 
accuracy after identifying the anomalies. 
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1. INTRODUCTION 

Hydrocarbon seepage systems are natural phenomena that occur when the hydrocarbon traps are not 
completely sealed. In such a case the hydrocarbons may escape from the subsurface source situated either 
offshore or onshore to the surface. It is estimated that about 75% of the world’s oil basins show hydrocarbon 
spills [1], [2]. In the onshore domain, the aspects of oil and gas seepage can be divided into two categories: 
a) macroseepage and b) microseepage. Macroseepages correspond to the accumulation of hydrocarbons 
(HCs) on the surface directly visible by the naked eye whereas microseepages are defined as traces of 
invisible light HCs in soils and which are only detectable by analytical methods. 

Hydrocarbon spills are a global problem, which can have a negative impact on human life and the 
environment. Therefore, early identification and taking necessary remedial measures are essential [3]. 
The hydrocarbons that appear on the surface cause oxidation and reduction reactions at sites or along migration 
paths. As a result, sediment and water anomalies occur [4]. In close-surface sediments, mineralogical changes 
and changes in electrical, magnetic and seismic property are observed [5]. Hydrocarbon pollution can affect 
vegetation and cause changes in vegetation (bio) physical parameters (pigments, water, and content), depending 
on the species, pollutant type, and exposition period [6]. 
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The hydrocarbon absorbs part of the solar radiation and emits it in the region 8 um and 14 um of the 
spectrum, as thermal energy. In infrared (IR) images on the basis of physics, every object with a temperature 
greater than absolute zero emits energy in the form of electromagnetic radiation [7], [8]. A black body is a 
hypothetical or model body that absorbs all the radiation falling on it and does not reflect or transmit it. 
It’s a perfect radiation emitter and a perfect absorber over all wavelengths. It takes values without dimensions 
ranging from (0 to 1) [9]. The thermal energy emitted from the black body is distributed spectrally depending 
on its temperature only. This radiation is referred to as thermal radiation. However, the amount of thermal 
radiation freed depends on the emitted surface of the body for any particular temperature and wavelength [10]. 

Remote sensing plays an important role in detecting oil spill response efforts to measure oil spills in 
large fields are often difficult because they require a very expensive and time-consuming processes such as 
inspecting large affected areas, affected area monitoring, and other field measurements. While remote 
sensing by satellite can provide information easily and effectively in large areas and at a lower cost [11]. 
One of the satellite sensors, landsat 8, was successfully launched in 2013 by National Aeronautics and Space 
Administration (NASA) and the United State Geological Survey (USGS); it is a ground-monitoring satellite 
that covers the same area once every 16 days. A scene of a landsat 8 image generally covers an area of 
185x185 km’. It is equipped with two sensors: the operational land imager (OLT) and the thermal infrared 
sensor (TIRS). OLI gathers data with a 30 m spatial resolution with eight bands situated in the visible, 
near-infrared, and shortwave infrared portions of the electromagnetic spectrum, as well as a 15 m 
panchromatic band. TIR radiance is measured at a spatial resolution of 100 m by TIRS, which employs two 
thermal bands situated in the atmospheric window between 10 and 12 um for continuous earth observations 
to provide for the estimation of land surface temperature (LST) [12], [13]. The landsat-8 satellite using both 
thermal image channels (i.e. 10 and 11) to provided temperature data. Band 10 was used to retrieve the 
Earth’s surface temperature using the algorithm described in the work steps. 

This study aimed to determine the quantities of hydrocarbons mixed with the soil as an indicator of 
its accumulation in the subsurface. Using a unique method for surveying and exploring an area of interest. 
This project aims to use remote sensing to detect the thickness of oil spills at different regions surrounding 
the oil fields. Various infrared bands will be examined to assess the accuracy of using them for determining 
the spill regions. And select the quantities of hydrocarbons mixed with the soil as an indicator of its 
accumulation in the subsurface. Using digital processing techniques for satellite visuals and Google images 
to deduce and identify surface evidence of oil spills. Detection of oil spills by remote sensing is a very 
important topic and many papers have been published in this field [14]. Such as radar images are able to 
operate at night or day and with clouds or fog. They can be used to detect oil spills in water, but they are 
limited to some monitoring conditions and have many limitations [15]. Also, the laser fluorosenseor is the 
most commonly utilized device. Due to its unique capacity to recognize oil on backdrops which include soil, 
water, ice and snow. The drawbacks are big size, weight and expensive cost [16]. 


2. RESEARCH METHOD 

The research area in this case study is located in southern Iraq between longitudes (48°30°E-46° 40'E) 
and latitude (31° 20° N-29° 50° N). This governorate is considered one of the important governorates in Iraq, 
as it is rich in oil fields. These include the Rumaila field, the Shaiba fields, the West Qurna field, and the 
Majnoon fields. This region has a high concentration of hydrocarbon (HC), particularly as an oil and gas 
producer.The occurrence of shale gas and HC seepage in this region indicates that HC is very potent. 
The study area is shown in Figure 1. 


2.1. Materials and methods 

In this study, landsat-8 data were used to estimate the land surface temperature based on thermal 
band 10 to calculate the brightness temperature and infrared band 5 and red band 4 to calculate normalized 
difference vegetation index (NDVI) of the study area. Knowing the NDVI value is essential to determine 
proportional vegetation and emissivity. The land surface temperature can be derived from these bands. 


2.2. The used data 

Landsat 8 is a United States ground-monitoring satellite (USGS) that covers the same area once 
every 16 days. Its data can be obtained free on the website of the USGS. Satellite imagery data were used for 
Al-Basrah region in October of 2020. The image was already rectified to WGS-1984-UTM-Zone_38N. 


2.3. The software’s used 


— ArcGIS 10.7: it is an advanced integrated system in geographic information systems, and it is produced 
by the American company environmental systems research institute (ESRI). ArcMap enables users to 
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view, modify, produce, and analyze geospatial data. This program is used for the final output of satellite 
images. 

— Surfer 16: surfer is a contouring and 3D surface mapping software program that runs under Microsoft 
Windows. The Surfer software quickly and easily converts data into outstanding contour, surface, 
wireframe, vector, image, shaded relief, and post maps. This program is used for the Interpolation map. 

— Matlab R2019b: Matlab allows matrix manipulations, plotting of functions and data, implementation of 
algorithms, creation of user interfaces, and interfacing with programs written in other languages. 
This program was used to solve a gradient equation and draw a map of it. 


Al-Basra 


47 0'0" 43°0'°0" 


35°0'0"N 


30°0'0"N 





40°0'0"E 45°0'0"E 


Figure 1. The study area located in South of Iraq, it is one of oil and gas onshore block in Iraq 


2.4. Methodology 
The method used in the proposed study to evaluate LST is depicted in Figure 2 and it is used to 
process landsat-8 data that were geometrically corrected. The proposed action steps are explained in detail in 
the following literature. 
—  Step-1: the proposed work’s first step is to use the following equation to convert digital number (DN) 
values to top of atmosphere radiance (TOA) for band 10 [12], [17]-[19]. 


Lı = M, X Qcal + A, — O; (1) 


Where, Qeai = quantized and calibrated standard product pixel values; L} = Top of atmosphere radiance 
(Watts / (m2 x sr x um));, A; = radiance add band 10; M, = radiance multiplicative band 10; O; = correction 
value for band 10 


—  Step-2: conversion TOA radiance to brightness temperate using thermal constant values in a metadata 
file and the (2). 


= ——_ — 273.15 (2) 
+1 


Where, BT = top of atmosphere brightness temperature (°C); K; and K2 are thermal band 10 constants that 

can be found in the satellite image’s metadata file. To convert the temperature from Kelvin to Celsius, the 

limit must be added (-273.15). 

—  Step-3: the NDVI is critical for identifying the study area’s various land cover types. The NDVI scale 
runs from -1 to +1. Through The difference between infrared band 5 and red band 4 can be calculated 
the NDVI of the images using (3). 


(NIR—RED) 


is (NIR+RED) (3) 


Where, NDVI = normalized difference vegetation index; NIR = near infra-red band; R = red band. 
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—  Step-4: proportional vegetation refers to estimating the proportion of an area It distinguishes vegetation 
from land cover by pixels. Pv can be calculated using (4). 


2 

NDVI-NDVI 

P, = (eos) (4) 
NDVIy—NDVI¢s 


Where, NDVIs is the NDVI of soil, NDVIv is the NDVI of vegetation. 

—  Step-5: the land surface emissivity (LSE) is a proportionality factor that measures the black body 
radiance to quantify emitted radiance and it is necessary to estimate LST. Furthermore, the LSE is 
highly affected by surface roughness, vegetation cover, and other factors: 


E = 0.004 x Pv + 0.986 (5) 


Where, E = the land surface emissivity and Pv = proportion of vegetation. 
—  Step-6: finally, the LST can be derived depending on the proportional vegetation and emissivity. 
Using the (6) [20]. 


BT 


ker (cia eer (6) 


tip) meal} 


Where, 7; = LST in Celsius (° C); 

A = the wavelength of band 10; 

£, = the emissivity; p = 1.438x107 mK 

—  Step-7 (interpolation): Kriging interpolation is based on statistical models (i.e. including auto-correlation) 
Geostatistical techniques can produce a forecast surface and also provide accurate forecasts. Kriging 
supposes that the distance or direction between the measured points reflects spatial correlation to 
explain the surface variation. The Kriging tool uses a mathematical relationship to determine the 
resulting values for each site within a given radius [21]. It is a multi-stage process that includes data 
exploration, analysis, surface creation, variogram modeling, and variance exploration. Kriging is best 
when you know that the data are spatially correlated or directed. Often it is used in the geology and soil 
sciences. Kriging is weighs its surrounding values to make a prediction for a non-rated location. 
Its general equation [22]. 


Z (So) = i AiZ (si) (7) 


Where: Z (s;) = the measured value at the i" location; 2; = unknown weight for the measured value at the i” 

location; So = the prediction location; N = the number of measured values. 

The weight in the ordinary Kriging process depends on the distance between the measured points 
and the expected location and the spatial relationship between the measured points around the location of the 
prediction, the parts that follow describe how to make a map of the prediction surface and a map of 
prediction accuracy using the general Kriging formula. To create a prediction surface map using the Kriging 
method there are two primary functions: 

a) Discover the rules of dependence 

b) Create the forecasts 

In order to perform these two functions, it is a two-stage process: 

a) It produces the covariance and variogram functions to assess the statistical dependency values that are 
dependent on the sample autocorrelation 

b) It forecasts unknown values [23] 

Summary the Kriging method uses the data twice: 

a) First, an estimate of the spatial relationship between the measured points 

b) Second, forecasting unknown values [24] 

—  Step-8 (gradient): the Image gradient is defined as the directional slope of an image’s tonal density and 
is used frequently for multivariate functions [25]. That incorporates pixel color changes on both the 
X-axes and y-axes. It is a vector of partial derivatives of all the variables [26]. Supposing that f (x, y) 
represents the color of the pixel at (x, y), the gradient vector of this location. It is defined as (8). 


| _ ee + ly) — f(x —1,y) 
f 


ECAN 
ED al loro esto- reyo (8) 
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The partial derivative of an x-axis is represented by the symbol (0f/ox). Measured as 
f (x+Ly) —f (x-J,y), similarly, the partial derivative of a y-axis is represented by the symbol (O//oy), measured 
as f (x,y+1) —f (%y-l). 
There are two important aspects to image gradient [27]: 
a) The magnitude is g = 49% + g3 
b) The direction is 0 = arctan (dy / Gx) 

The gradient of an image is a measurement of how it changes. It provides two types of information. 
The gradient magnitude tells us how fast the image is changing. While the gradient direction tells us the 
direction in which the image is changing [28]. This vector’s length supplies the gradient value. While it gives 
its direction the gradient direction. Since the gradient can vary at each image location, we use a different 


vector to represent it [29]. 
| Landsat $ | 
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Figure 2. Flow chart diagram 


3. RESULTS AND ANALYSIS 

The findings revealed that there are 60 thermal anomaly locations in the study area and that the lowest 
and highest temperature value is 24.2 °C to 91.11 °C. Three different sites were chosen in terms of shape, area, 
and effect. Which is located in the northern and southern Rumaila and Zubair fields, as in Figure 3, Figure 4 and 
Figure 5. 

After identifying the anomalous areas. The spatial distribution of the temperature in the region was 
obtained by applying the Kriging method (which is one of the geo-statistical interpolation methods. For this 
reason, interpolation methods are used to getting complete surface information). Spatial interpolation is the 
process of calculating an unknown value using a group of sample points with known values dispersed 
throughout a region. The distance between the cell with uncertain value and the sample cells adds to the 
estimate of its ultimate value. The objective of spatial interpolation is to generate a surface that best 
represents the sampled event. The first site’s interpretation of interpolation techniques was displayed. 
The shape is semi-oval, with the darker blue color in the center representing the maximum temperature value, 
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and the further away from the middle the color tends to light because the temperature gradually decreases. 
It is one source of hydrocarbon seepage, and it represents the well, which appears in dark blue, and the light 
blue color indicates the spread of hydrocarbons around it. As in Figure 6, while at the second site there are 
two hydrocarbon seepage close together, one strong appears dark blue and the other is weak, as in Figure 7. 
Finally, the third site shows more than one source of hydrocarbon seepage spreading over a wide distance. 
It also appears in Figure 8 in dark blue. A gradient determination method has been applied to get more 
information about its location, extension, size, shape, and trends. It also appears in Figure 9, Figure 10, and 
Figure 11 where the direction of the arrow indicates areas of leakage. The result of gradient and interpolation 
methods was effective to detect evidence of hydrocarbon estimate significant areas that have hydrocarbon 
seepage in the study area which may be promising for exploration and production of oil and gas. 
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Figure 3. Image Google Earth view location of hydrocarbon seepage for Southern Rumaila oil field 
in the anomalous area-1 
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Figure 4. Image Google Earth view location of hydrocarbon seepage for Zubair oil field in the anomalous 
area-2 
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Figure 5. Image Google Earth view location of hydrocarbon seepage for Northern Rumaila oil field in the 
anomalous area-3 
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Figure 6. Interpolation of area-1 using Kriging Figure 7. Interpolation of the H.C distribution for 
method area-2 using Kriging method 
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Figure 8. Interpolation of the H.C distribution for Figure 9. The gradient map for Southern Rumaila oil 
area-3 using Kriging method 
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Figure 10. The gradient map for Zubiar oil field Figure 11. The gradient map for Northern Rumaila 
oil field 


4. CONCLUSION 

This paper presents a new method for determining the locations of hydrocarbon leaks on the surface 
using the temperature data obtained from the landsat-8 satellite images to determine the extent of the oil spill 
statistically by using spatial interpolation and gradient techniques. This study illustrates that landsat-8 TIR 
can be used for the exploration of gas and oil resources by HC micro seepage detection in the region. 
The detection of gas and oil resources is conducted by temperature data. The findings of this study show that 
the HC micro-seepage range spreads in the region around the oil field. The results showed high accuracy in 
determining the locations of oil spills and the size of them after identifying the anomalies areas The results 
achieved were very rational and showed high accuracy in determining the source of the leak when compared 
to Google Earth images, and we note the presence of a flame at the site of the leak. The proposed approach in 
this study would help determine the shape of the anomaly and reveal the spatial relationships between the 
mapped anomaly and possible subsurface. 
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